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The electrical contact between an embed- 
ded microelectrode and a cultured neu- 
ron depends on the geometry of the 
neuron-electrode interlace. The contact is 
improved when the electrode is covered, or 
sealed, completely by the neuron. In this 
article, the finite element method is pro- 
posed as a tool for modeling the electrical 
properties of the neuron-electrode inter- 
face. This method permits numerical solu- 
tions of volume conductor problems for a 
\ ariety of geometries, without prior restric- 
tion of the current paths. Simulations are 
focused on the influence of the geometry 
on the transfer of an extracellularly applied 
stimulus current to the neuron and on the 
sealing resistance. A comparison is also 
made between finite element modeling and 
lumped circuit modeling. 

Overview 

One of the objectives of the 
neurotechnolgy program in our group is to 
establish a neuro-electronic interface with 
the human nervous system on the cellular 
level [1]. For the development of this 
implantable device, a two-step strategy is 
followed. In the first step, a large number 
of electrical contacts arc made between an 
array of microelectrodes and a layer of 
cultured neurons. In the second step, after 
implantation, synaptical contacts will be 
formed between the neuronal culture and 
the nervous system. 

The electrical contact between 
microelectrodes and cultured neurons can 
be studied in detail using planar substrate 
multielectrode arrays (MEAs). In Fig. 1(a), 
one of the MEAs that is currently used by 
our group is depicted. It consists of a hexag- 
onal array of 61 microelectrodes, each with 
a diameter of 10 urn. In this example, the ar- 
ray is covered with neuronal cells that are 
isolated from dorsal root ganglions of neo- 
natal rats. In other groups, similar MEAs are 
widely used with increasing success for 
long-term recording and stimulation of neu- 
ral networks in culture [2-4]. From the large 




number of cell-size electrodes, each con- 
tacting one or a few individual cultured 
neurons, a diversity of shapes and a wide 
range of amplitudes of extracellularly re- 
corded action potentials have been re- 
ported [5, 6]. Hence, a variability exists in 
the behavior of the electrical contacts be- 
tween the neurons and microelectrodes, 
which also has consequences for the stimu- 
lation of individual neurons. 

To establish more reliable and 
well-defined neuron-electrode contacts, a 
proper understanding is required of the 
way electrical signals are transferred be- 
tween a neuron and an electrode. Ba- 
sically, this transfer is mediated by the 
volume conducting properties of the 
extracellular space — the culture medium. 
This volume conductor and its electrical 
interactions with the electrode and the 
neuron is the primary subject of this 
article: it is called the neuron-electrode in- 



terface (not to be confused with the 
neuro-electronic interface). Two exam- 
ples of neuron-electrode interfaces are de- 
picted in Figs. 1(b) and 1(c). Although 
electrodes are often covered only partially 
with neurons, sometimes the neurons 
cover the electrode completely. 

When a neuron covers, or seals, an elec- 
trode completely, a specific part of the in- 
terface becomes important for the 
neuron-electrode contact: the sealing gap 
between the neuronal membrane and the 
substrate (Fig. 2). Experimental results 
show that in this case the amplitude of re- 
corded action potentials increases and the 
stimulation threshold decreases [6, 7]. 
These improvements of the contact are as- 
cribed to the decreased resistance of the 
sealing gap. Hence, the electrical proper- 
ties of the neuron-electrode interface are 
related to its geometry; i.e., size, shape and 
relative positions of neuron andelectrode. 

For interface geometries of complete 
sealing, it has become common practice to 
model the neuron-electrode contact with 
simple lumped circuits [5-7], based on the 
two major current pathways that exist in 
these cases: the leakage current and the 
current through the cell (Fig. 2). The pa- 
rameters of these circuits are approxi- 
mated from the interface geometry. 
However, the validity of this method is 
questionable for at least two reasons. 
First, major simplifications arc made in 
the representation of the electrical proper- 
ties of the interface by the definition of 
only two current paths. Second, it is diffi- 
cult to compute the component values 
from the geometry, unless the represented 
geometry is very simple (e.g., 
axisymmetrical) or can be approximated 
by discretization into simple elements. 

Methods 

Parametrical Geometry 

Although the finite element method 
permits modeling of almost any geome- 
try, a parametrical geometry is defined to 
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gain insight into the relationship between 
geometry of the neuron-electrode inter- 
face and its electrical behavior. This 
parametrical geometry is depicted in Fig. 
3. The neuron is modeled as a circular 
soma of radius r c , with a parabolic height 
profile h c . The axon and dendrites of the 
neuron are not included in this version of 
the model, since we are primarily inter- 
ested in situations in which neurons cover 
electrodes completely or partially. In each 
simulation, the radius of the neuron is var- 
ied from r c = 7 tor c = 20 jam in four steps. 
These values are representative of a wide 
variety of vertebral neuronal types, such 
as cortical neurons, dorsal root ganglion 
neurons, and spinal cord neurons [8-10]. 
The height of the neuron is taken as h c = 5 
Jim, based on the range of neuronal diam- 
eters and the assumption of a certain flat- 
tening of the neuron in culture, due to its 
adhesion to the substrate. The neuron is 
positioned with an eccentricity x c on top 
of an electrode with radius r e = 5 jam, 
which is located in the bottom of a well 
with a depth of h e = 0.8 (im. A sealing gap 
of thickness d g = 10 nm is modeled be- 
tween the soma and the substrate, repre- 
senting a strong neuron-substrate 
adhesion [11]. 

Finite Element Model 

The geometry of the neuron-electrode 
interface is implemented in ANSYS Fi- 
nite Element Software (SAS IP Inc. Hous- 
ton, TX) release 5.4 [Fig. 4(a)]. The 
medium surrounding the interface is rep- 
resented by a volume conductor of con- 
ductivity ct =1.65 S/m. The 
intracellular volume is modeled with a 
conductivity of o = 1.43 S/m. The 
thickness of the sealing gap </,, is two or- 
ders of magnitude smaller than the other 
geometry parameters. Implementing the 
sealing gap with its actual thickness 
would require a mesh with either imprac- 
tical element shapes, or an unfeasibly 
large number of elements. Therefore, the 
sealing gap is implemented with a fixed 
thickness of d = 0.5 (tm, and its 

conductivity is compensated for the dif- 
ference between modeled and imple- 
mented thickness, according to: 




All volumes in the model are meshed 
and filled with tetrahedral-shaped volume 
elements (SOLID98), which permit nu- 
merical solutions of the Poisson equation: 
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V-(ctW) = 0 (2 ) 

with V the electrical potential and a the 
conductivity of the culture medium, the 
intracellular fluid, or the fluid in the seal- 



ing gap. The size of the elements varied 
from 0.5 um in the sealing gap to several 
microns far away from the neuron-elec- 
trode interface. 





The electrode 
impedance has no 
effect on the stimulus 
transfer and therefore 
the electrode 
impedance is not 
required in this model. 



The electrical properties of the 
neuronal membrane are assumed to be 
highly anisotropic. Due to the ion chan- 
nels, current can only flow in the direction 
perpendicular to the surface of the mem- 
brane. Hence, there is no need for model- 
ing the membrane using volume elements, 
and the membrane can be represented by a 
set of extracellular ant! intracellular mem- 
brane nodes [Fig. 4(b)]. Each intracellular 
membrane node corresponds to one 
extracellular membrane node and is con- 
nected to it by a resistive circuit element. 
The resistance of this element is com- 



puted using the membrane conductivity 
and the membrane area represented by the 
connected node. Membrane capacity is 
left out of the model, since pilot simula- 
tions with a membrane capacity of 1 
p.F/cm 2 revealed that the time constant of 
the membrane potential is less than a 
microsec (results not presented). The 
membrane potential is computed from the 
potential difference between an intra- 
cellular membrane node and its corre- 
sponding extracellular membrane node. 
No active (e.g., Hodgkin-Huxley) proper- 
ties of the membrane are modeled. 

Since current stimulation is used, the 
electrode impedance has no effect on the 
stimulus transfer and therefore the elec- 
trode impedance is not required in this 
model. All nodes on the surface of the elec- 
trode are coupled in the voltage domain, 
and a stimulus current is appl ied to the total 
electrode surface. The nodes at the outer 
boundary of the modeled part of the culture 
medium are set to zero potential, represent- 
ing a counter electrode that is far away 
from the neuron-electrode interface. 

In each simulation, both the stimulus 
transfer and the resistance of the sealing 
gap are computed. The stimulus transfer 
is defined as the maximum depolariza- 
tion of the upper membrane, due to a 
stimulus current of 7 = 1 nA through 
the electrode. When the neuron com- 
pletely covers the electrode, this stimulus 
current will depolarize the upper mem- 
brane by several millivolts. This permits 
easy interpretation of the results with re- 



spect to the stimulus threshold for 
eliciting an action potential. 

From previous work, it is known that at 
complete sealing, the current through the 
cell is negligible compared to the current 
through the sealing gap. The potential at 
the electrode surface is therefore mainly 
determined by the sealing resistance, 
which is thereby computed as: 



7 stim 



(3) 



When the sealing becomes defective, a 
sealing resistance no longer exists. Instead, 
the value computed from Eq. (3) must be 
interpreted as a spreading resistance. 

Results 

Complete and Defect Sealing 

In practice, cultured neurons rarely 
cover an electrode completely. Therefore, 
the effects of varying eccentricity of a neu- 
ron with a radius of;;. = 7 urn are studied 
using eccentricities of x c = 0 (complete 
sealing) and x c = 9 um (defective sealing). 
In Fig. 5, the computed potential distribu- 
tion due to a 1 nA stimulus current is de- 
picted in a cross section of the model. The 
equipotential lines in the medium are 
marked by the boundaries between two ad- 
jacent colored regions. At complete seal- 
ing | Fig. |, these equipotential lines are 
concentrated in the sealing gap, indicating 
a potential drop over the sealing gap. A 
maximum depolarization of 1.02 mV 
(stimulus transfer) occurs at the upper 
membrane, and a hyperpolarixationof 2.25 
mV occurs at the lower membrane. When 




-- Membrane Node 



4. (a) 3-D visualization of the neuron electrode interface geometry, as implemented in the ANSYS finite element software, (b) 
The membrane is represented by a set of resistive circuit elements, connecting the extracellular and intracellular membrane 
nodes. 
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the sealing is defective [Fig. 5(b)], the 
equipotential lines are spread up more 
widely, indicating a leakage current into 
the medium. The changes in membrane po- 
tential are reduced to the microvolt range. 
From the electrode potential V e , and the 
stimulus current at x c = 0, a sealing resis- 
tance of 3.28 MQ is computed and at 
x c = 9 Lim a spreading resistance of only 
36.6 k£2 is computed. These results reveal 
that sealing resistance is strongly related to 
the stimulus transfer. The transition from 
complete to defective sealing causes a 
drastic decrease in both the stimulus trans- 
fer and the sealing resistance. 

Neuron-Substrate Adhesion 

For culturing monolayers of neurons, a 
proper adhesion of the neurons to the sub- 
stratum is desired. Often, adhesion pro- 
moters are used for this purpose, such as 
poly-l-lysin, laminin, and polye- 
thyleneimine. Due to a stronger adhesion 
to the substrate, neurons generally will 
spread out and become flatter. Hence, ad- 
hesion will affect the radius and the height 
of the neuron. Furthermore, the electrical 
properties of the sealing gap can be af- 
fected. These properties are the subject of 
some consideration in the literature. Gen- 
erally, a bulk resistivity is assumed in the 
gap and the thickness of the gap is assumed 
to be very small, which results in sealing 
resistances of several megohms. 

In Fig. 6, the simulated effects of varia- 
tions in neuronal radius and height, and 
also of the thickness of the sealing gap, are 
presented for cases of complete sealing 
(x c = 0). The stimulus transfer [Fig. 6(a)] 
and the sealing resistance [Fig. 6(b]) are 
plotted versus the radius of the neuron. 
The solid line represents the results from 
simulation with normal height and thick- 
ness of the sealing gap. Both stimulus 
transfer and sealing resistance increase 
with increasing radius of the neuron. The 
dashed lines represent results from simu- 
lations in which one geometry parameter 
(h c ord g )is varied. 

Although the height of the neuron is 
varied extensively to h c = 1 um and to 
h c = 9 urn, it does not have a large effect on 
the stimulus transfer, and its effect on the 
sealing resistance is even smaller. How- 
ever, when the thickness of the sealing gap 
is doubled from d g = lOnmto^ = 20 nm, 
both the stimulus transfer and the sealing 
resistance are reduced by a factor of two. 
The same effect appears when the sealing 
gap is halved from d g = 10 nm to d g = 5 
nm: both the sti i nulus transfer and the seal- 



ing resistance are doubled, regardless of 
the size of the neuron. 

Finite Element Model 
versus Circuit Model 

As stated at the start, finite element 
modeling is proposed in this article as an 
alternative for lumped circuit modeling. 
In order to make a comparison between 




these two modeling methods, the 
parametrical geometry of the neu- 
ron-electrode interface is also used for 
computation of the components of a 
lumped circuit. Just as in the finite ele- 
ment model, only resistive membrane 
properties are modeled. Furthermore, 
since current stimulation is used, the elec- 
trode impedance is left out of the model. 
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6. Effect of variation of several model parameters on upper membrane depolarization 
and on sealing resistance. The results from simulations with the initial parameter val- 
ues (h c = 5 iim,d g = 10 nm) are plotted with a solid line. The dashed lines represent 
results from simulations in which one of the parameters is varied, (a) Effect of varia- 
tion of the height of the cell ( V h c =1 um, □ h c = 9 um) and the thickness of the seal- 
ing gap (o d g =5 nm, x d g = 20 nm) on the upper membrane depolarization, (b) 
Effect of variation of the height of the cell ( Vh c = 1 um, □ h c = 9 nm) and the thick- 
ness of the sealing gap (o d g =5 nm, x d g =20 nm) on the sealing resistance. 
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5. Computed potentials due to a 1 nA stimulus current through the microelectrode 
in a cross section of the model. The equipotential lines are marked by the bound- 
aries between two colored regions. 
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The resulting circuit model consists of 
three resistive components: the sealing re- 
sistance and the upper and lower mem- 
brane resistances (Fig. 7). 

The depolarization of the upper mem- 
brane due to a stimulus current through 
the electrode can be written as: 

y up = * seal 7sdm 

^seal + ^low + ^up (4) 

When the neuron is centered on the 
electrode (x c = 0), the geometry parame- 
ters can be used for computing the circuit 
components: 



K-°mem- r l (6) 



"up T 

K-Gmem-fiKyrc (7) 

with the radius of the lower membrane 
r t = r e and f(h c ) = 1.43, a correction fac- 
tor for the area of the upper membrane, 
due to the height of the neuron. 

An interesting test case for the compar- 
ison between finite element modeling and 
lumped circuit modeling is variation of 
the radius of the electrode. The electrode 
radius is a design parameter of the MEA 
and can be adjusted for optimization of the 
neuron-electrode contact. A decrease of 
the radius of the electrode r e will have two 
contradictory effects on the stimulus 
transfer: On one hand, the length of the 
current path through the sealing gap is in- 
creased, thereby increasing the sealing re- 
sistance and the stimulus transfer. On the 
other hand, the area of the lower mem- 



brane (covering the electrode surface) is 
decreased, thereby decreasing the current 
through the cell and, hence, the stimulus 
transfer. Depending on which of these two 
effects is the strongest, the stimulus trans- 
fer may increase or decrease with decreas- 
ing electrode radius. Also, an optimal 
electrode radius may exist with respect to 
the stimulus transfer. Obviously, the net 
effect of a decreasing electrode radius on 
the stimulus transfer must be obtained by 
proper modeling of the neuron-electrode 
interface, which makes this an ideal test 
case for a comparison between the two 
modeling methods. 

In Fig. 8, the stimulus transfer and 
sealing resistance, as computed with both 
the finite element model and the circuit 
model, are plotted versus the radius of the 
neuron. In these simulations, the electrode 
radius is decreased in three steps from 
r e = 5 Jim to r e = 3 urn. For the stimulus 
transfer, there appears to be a striking dif- 
ference between the two models [(Fig. 
8(a)]. The results of the circuit model sug- 
gest that an optimal stimulus transfer ex- 
ists for a certain radius of the neuron, and 
that this optimum depends on the size of 
the electrode. The results from the finite 
element model, however, do not show this 
optimum. Instead, the finite element 
model predicts a higher stimulus transfer. 

The results in Fig. 8(b) show that this 
disagreement cannot be explained by dif- 
ferences in the computed sealing resis- 
tance. Instead, it arises from the membrane 
over the sealing gap, which is not repre- 
sented in the circuit model. Due to the po- 
tential drop that exists over this sealing 
membrane, a radially decreasing mem- 
brane current arises from the sealing gap 
into the intracellular space (Fig. 9). This 
current is negligible compared to the leak- 
age current in the sealing gap and will not 
affect the computation of the sealing resis- 
tance. However, it is not negligible com- 
pared to the current through the part of the 
lower membrane that covers the electrode. 
Hence, this membrane current contributes 
significantly to the stimulus transfer. 

Nonparametrical Geometry 

The finite element method permits 
modeling of almost any geometry and is 
not restricted to parametrical geometries, 
as defined above. Therefore, more realis- 
tic shapes of neurons can be modeled eas- 
ily, which is very advantageous when it 
comes to experimental characterization of 
the electrical properties of the neu- 
ron-electrode interface. As an illustration, 
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7. Circuit model of the neuron-electrode interface, used for comparison with the fi- 
nite element model. The electrode impedance is left out of the model and only resis- 
tive membrane properties are modeled. All component values are computed from 
the geometry parameters (Fig. 3). 




8. Comparison between simulation results from circuit mode! (solid lines) and finite 
element model (dashed lines), (a) Upper membrane depolarization, plotted versus 
the radius of the neuron while the radius of the electrode is varied (+ r e =5 um, o 
r e = 4 um, Vr e = 3 urn) (b) Sealing resistance, plotted versus the radius of the neuron 
while the radius of the electrode is varied (+ r e =5 um, o r e = 4 urn, Vr c = 3 um). 
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the two examples of neuron-electrode in- 
terfaces from the beginning of this article 
[Figs. 1(b), Fig. 1(d)] are used for recon- 
struction of the interface geometry. The 
circumferences of neuronal soma and 
electrode are marked on a digitized micro- 
photograph. From these contours, the ge- 
ometries are reconstructed [Fig. 10(a)]. 
The assumption of a parabolic height pro- 
file is applied again, since the above pre- 
sented simulations revealed that the 
influence of height was not very large. 
The maximum height was taken as a quar- 
ter of the maximal diarnetei'of the neuron. 

Based on the reconstructed geometry, 
mesh generation is performed in ANSYS. 
Again, a stimulus current is applied 
through the electrode, and the membrane 
potentials are computed. In Fig. 10(b), the 
membrane potential distribution is de- 



picted on the model that is superposed on 
an image of the electrode. A view, from 
below, in Fig. 10(c), reveals that most of 
the membrane polarization is localized at 
the site of electrode and in the sealing gap. 
The computed potentials can be related 
directly to experimentally derived poten- 
tials of the electrode and the intracellular 
space, which permit identification of the 
electrical properties of the neu- 
ron-electrode interface. 

Conclusion 

The finite element model offers sev- 
eral advantages above lumped circuit 
modeling. It permits linking of the electri- 
cal properties of the neuron-electrode in- 
terface to a variety of more or less 
complex geometries, including cases of 
both complete and defective sealing. All 




currents and potentials in the neu- 
ron-electrode interface are computed in- 
herently by generating a numerical 
solution of the finite element problem, 
which obviates the need for simplifying 
the electrical representation of the geome- 
try. Furthermore, this model can be easily 
extended with the implementation of ac- 
tive membrane properties and axons, 
which allow simulation of extracellular 
stimulation as well as recording. 

The applications of this modeling 
method are numerous. For example, using 
geometrical data it can be used for electri- 
cal characterization of the sealing gap and 
the effects of several adhesion promoters 
on sealing resistance. It also permits mod- 
eling of the influence of local (active) 
membrane properties, such as ion-specific 
channel densities, on extracellular record- 
ing and stimulation. Finally, it would 
greatly facilitate the optimization of the 
MEA substrate and electrode geometries 
for improved neuron-electrode contact. 
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